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What is claimed is: 

1. A method of estimating physical parameters of porous rocks of a subterranean formation 
containing a compressible reservoir fluid comprising the steps of: 

(a) injecting an injection fluid into the subterranean formation at an injection pressure 
exceeding the subterranean formation fracture pressure; 

(b) shutting in the subterranean formation; 

(c) gathering pressure measurement data over time from the subterranean formation during 
shut-in; 

(d) transforming the pressure measurement data into corresponding adjusted pseudopressure 
data to minimize error associated with pressure-dependent reservoir fluid properties; and 

(e) determining the physical parameters of the subterranean formation from the adjusted 
pseudopressure data. 

2. The method of claim 1 wherein a plot of the adjusted pseudopressure data over time is a 
straight line with a slope m M and an intercept b M , wherein m M is a function of permeability k 
and b M is a function of fracture-face resistance R 0 . 

3. The method of claim 2 wherein the adjusted pseudopressure data used in the transforming 
step are derived using following equation: 
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pdp_ ^ w herein 
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Jo Mg c t 


M 


= average viscosity, m/Lt, cp 




= gas viscosity, m/Lt, cp 


p 


= pressure, m/Lt 2 , psi 


p 


= average pressure, m/Lt 2 , psi 


Pa 


= adjusted pseudopressure variable, m/Lt 2 , psi 


Pw 


= wellbore pressure, m/Lt 2 , psi 


PLfD 


= dimensionless pressure in a hydraulically fractured well, dimensionless 


c t 


= total compressibility, Lt /m, psi" 




= average total compressibility, Lt 2 /m, psi" 1 
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4. The method of claim 3 wherein the straight line is defined by the equation: 
(y a )« = m M (x a )„ + b M , where 

/„ v _ (Pa)n-Par 
\sa )n ~ I — I , 



[Pa<J>)] jA -[Pa(P)]j 



tj-tj-\ 



, and 
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(^a)n V Vwe ) 
" Kd a )j-(d a )j A ] 



j=ne+3 
c a2 
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wherein 

d a 

Par 
Paw 
tn 
tne 

(y a )n 



a first before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , 

„ -1/2 \I2 

psi *cp 

a second before-closure pressure-transient analysis adjusted variable, 

J/t 2/7/2 -3/2 1/2 

m/Lt , psi *cp 

before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr 

adjusted pressure variable difference, m/Lt 2 , psi 

adjusted reservoir pressure variable, m/Lt 2 , psi 

wellbore adjusted pressure variable, m/Lt 2 , psi 

time at timestep w, t, hr 

time at the end of an injection, t, hr 

before-closure pressure-transient analysis adjusted variable, dimensionless 
before-closure pressure-transient analysis adjusted variable, dimensionless 



5. The method of claim 4 wherein the first and second before-closure pressure-transient 
analysis variables are defined as: 



< a \ s ; 



and 



5.615 0 -g 



24 
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wherein 

^ = porosity, dimensionless 

B g = gas formation volume factor, dimensionless, bbl/Mscf 

B g = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L t , psi/ft 

w L = fracture lost width, L, ft 

6. The method of claim 5 wherein the transforming step is iterated with a value of n varying 
from ne+l to a maximum value n max and for each couple of coordinates {(y a ) n , OOn} plot the 
graph (y a ) n versus (x a )n to determine the slope m M and the intercept b\{> 

wherein 

ne = number of measurements that corresponds to the end of an injection 
n ma x = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure 

7. The method of claim 6 wherein the permeability k and the fracture-face R 0 are determined by 
the following equations: 

k = 



(141.2)(2)(0.02878)(24) 1 



5.615 r pSf m M 



^ = U\TJ(24) r P S f tnebM 



8. The method of claim 6 wherein the permeability k and the fracture-face Ro are determined by 
the following equations: 

-i2 



(141.2)(2)(0.02878)(24) 1_ 

5.615 r p Sftn M 



wherein 

o) = natural fracture storativity ratio, dimensionless. 

9. The method of claim 1 wherein the injection fluid is a liquid, a gas or a combination thereof. 
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10. The method of claim 9 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation. 

11. The method of claim 1 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 
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12. A method of estimating physical parameters of porous rocks of a subterranean formation 
containing a compressible reservoir fluid comprising the steps of: 

(a) injecting an injection fluid into the subterranean formation at an injection pressure 
exceeding the subterranean formation fracture pressure; 

(b) shutting in the subterranean formation; 

(c) gathering pressure measurement data over time from the subterranean formation during 
shut-in; 

(d) transforming the pressure measurement data into corresponding adjusted pseudopressure 
data and time into adjusted pseudotime data to minimize error associated with pressure- 
dependent reservoir fluid properties; and 

(e) determining the physical parameters of the subterranean formation from the adjusted 
pseudopressure and adjusted pseudotime data. 

13. The method of claim 12 wherein a plot of the adjusted pseudopressure data over time is a 
straight line with a slope m M and an intercept b M , wherein m M is a function of permeability k 
and bfA is a function of fracture-face resistance Rq. 

14. The method of claim 13 wherein the adjusted pseudotime and adjusted pseudopressure data 
used in the transforming step are respectively determined by the following equations: 




77 c. f(Pw)n 
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, wherein 
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1* = average viscosity, m/Lt, cp 

fi g = gas viscosity, m/Lt, cp 

p = pressure, m/Lt 2 , psi 

p = average pressure, m/Lt 2 , psi 

p a = adjusted pseudopressure variable, m/Lt 2 , psi 

p w = wellbore pressure, m/Lt 2 , psi 

PL f D = dimensionless pressure in a hydraulically fractured well, dimensionless 

c t = total compressibility, Lt 2 /m, psi" 1 
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= average total compressibility, Lt 2 /m, psi' 



15. The method of claim 14 wherein the straight line is defined by the equation: 
(yap )n= b M + m M i x ap )n , where 

\ = (Pa)n ~ Par 
\yap)n - 
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wherein 

c apl= c al 

dap 

&Pa 
Par 
Paw 
*n 

(Qn 

( x ap)n 

(yap)n 



a first before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , 

•1/2 1/2 

psi -cp 

a second before-closure pressure-transient analysis adjusted variable, 

2 /T 2 f 7/2 -3/2 1/2 

m/Lt , psi cp 

before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr, 

with adjusted pseudotime variable 

adjusted pressure variable difference, m/Lt 2 , psi 

adjusted reservoir variable pressure, m/Lt 2 , psi 

wellbore adjusted pressure variable, m/Lt 2 , psi 

time at times tep w, t, hr 

time at the end of an injection, t, hr 

adjusted time at timestep w, t, hr 

before-closure pressure-transient analysis adjusted variable, dimensionless 
before-closure pressure-transient analysis adjusted variable, dimensionless. 



16. The method of claim 15 wherein the first and second before-closure pressure-transient 
analysis variables are defined as: 
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and 



5.615 Bg K . 

24 



wherein 

^ = porosity, dimensionless 

^ = gas formation volume factor, dimensionless, bbl/Mscf 

5^ = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L 2 t 2 , psi/ft 

w L - fracture lost width, L, ft. 

17. The method of claim 16 wherein the transforming step is iterated with a value of n varying 
from ne+\ to a maximum value n max and for each couple of coordinates {(y ap )n, (xa P )n} plot the 
graph {y ap ) n versus (x op ) n to determine the slope /has and the intercept , 

wherein 

ne = number of measurements that corresponds to the end of an injection 
n ma x = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure. 

18. The method of claim 17 wherein the permeability k and the fracture-face R 0 are determined 
by the following equations: 

l2 



k = 



(141 .2X2X0-02878)(24) 1 

5.615 r pS/m M 



^ =: ^k^) rpSftnebM - 



19. The method of claim 17 wherein the permeability k and the fracture-face R 0 are determined 
by the following equations: 

(141.2)(2)(0.02878)(24) 1 

0)k = 



(141.2)(2)(0.02878)(24) 



5.615 



r p S f m M 



5.615 



141.2^(24) 



r p S f l ne b M 
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wherein 

o) = natural fracture storativity ratio, dimensionless. 

20. The method of claim 12 wherein the injection fluid is a liquid, a gas or a combination 
thereof. 

21. The method of claim 20 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation. 

22. The method of claim 12 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 
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23. A method of estimating permeability k of porous rocks of a subterranean formation 
containing a compressible reservoir fluid comprising the steps of: 

(a) injecting an injection fluid into the subterranean formation at an injection pressure 
exceeding the subterranean formation fracture pressure; 

(b) shutting in the subterranean formation; 

(c) gathering pressure measurement data over time from the subterranean formation during 
shut-in; 

(d) transforming the pressure measurement data into corresponding adjusted pseudopressure 
data to minimize error associated with pressure-dependent reservoir fluid properties; and 

(e) determining the permeability k of the subterranean formation from the adjusted 
pseudopressure data. 

24. The method of claim 23 wherein a plot of the adjusted pseudopressure data over time is a 
straight line with a slope thm which is a function of permeability k. 

25. The method of claim 24 wherein the adjusted pseudopressure data used in the transforming 
step are derived using the following equation: 
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= average viscosity, m/Lt, cp 


fig 


= gas viscosity, m/Lt, cp 


p 


= pressure, m/Lt 2 , psi 


p 


= average pressure, m/Lt 2 , psi 


Pa 


= adjusted pseudopressure variable, m/Lt 2 , psi 


Pw 


= wellbore pressure, m/Lt 2 , psi 


PLfD 


= dimensionless pressure in a hydraulically fractured well, dimensionless 


c t 


= total compressibility, Lt 2 /m, psi" 1 




= average total compressibility, Lt 2 /m, psi" 1 . 



26. The method of claim 25 wherein the straight line is defined by the equation: 
(y a )n = m M ( x a )n + b M , where 
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i \ = (Pg)n Par 
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wherein 

APa 

Par 

Paw 

t n 

*ne 

( x a)n 

(y a ) n 



a first before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , 
psi ,/2 cp ,/2 

a second before-closure pressure-transient analysis adjusted variable, 

2 /T 2/7/2 -3/2 1/2 

m/Lt , psi xp 

before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr 

adjusted pressure variable difference, m/Lt 2 , psi 

adjusted reservoir pressure variable, m/Lt 2 , psi 

wellbore adjusted pressure variable, m/Lt 2 , psi 

time at timestep w, t, hr 

time at the end of an injection, t, hr 

before-closure pressure-transient analysis adjusted variable, dimensionless 
before-closure pressure-transient analysis adjusted variable, dimensionless. 



27. The method of claim 26 wherein the first and second before-closure pressure-transient 
analysis variables are defined as: 



. . = IEL. 



and 



5.615 _ B g ui 
wherein 

<p = porosity, dimensionless 
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B g = gas formation volume factor, dimensionless, bbl/Mscf 

B g = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L 2 t 2 , psi/ft 

w L = fracture lost width, L, ft. 

28. The method of claim 27 wherein the transforming step is iterated with a value of n varying 
from ne+\ to a maximum value n max and for each couple of coordinates {{y a ) n , (x a )n} plot the 
graph (y a ) n versus (x a ) n to determine the slope yhm, 

wherein 

ne = number of measurements that corresponds to the end of an injection 
n ma x = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure. 

29. The method of claim 28 wherein the permeability k is determined by the following 
equation: 

i2 



k = 



(141.2)(2)(0.02878)(24) 1 

5.615 r pSf m M 



30. The method of claim 28 wherein the permeability k is determined by the following 
equation: 

hai ?v?vn m%7%vjd\ i 

cok = 



(141.2)(2)(0.02878)(24) 1 

5.615 r pSf m M 



wherein 

co = natural fracture storativity ratio, dimensionless. 

31. The method of claim 23 wherein the injection fluid is a liquid, a gas or a combination 
thereof 

32. The method of claim 31 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation. 
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33. The method of claim 23 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 



L 
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34. A method of estimating permeability k of porous rocks of a subterranean formation 
containing a compressible reservoir fluid comprising the steps of: 

(a) injecting an injection fluid into the subterranean formation at an injection pressure 
exceeding the subterranean formation fracture pressure; 

(b) shutting in the subterranean formation; 

(c) gathering pressure measurement data over time from the subterranean formation during 
shut-in; 

(d) transforming the pressure measurement data into corresponding adjusted pseudopressure 
data and time into adjusted pseudotime data to minimize error associated with pressure- 
dependent reservoir fluid properties; and 

(e) determining the permeability k of the subterranean formation from the adjusted 
pseudopressure and adjusted pseudotime data. 

35. The method of claim 34 wherein a plot of the adjusted pseudopressure data over adjusted 
pseudotime data is a straight line with a slope m M which is a function of permeability k. 



36. The method of claim 35 wherein the adjusted pseudotime and adjusted pseudopressure data 
used in the transforming step are respectively determined by the following equations: 
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= average viscosity, m/Lt, cp 




= gas viscosity, m/Lt, cp 
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= pressure, m/Lt 2 , psi 
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= average pressure, m/Lt 2 , psi 


Pa 


= adjusted pseudopressure variable, m/Lt 2 , psi 


Pw 


= wellbore pressure, m/Lt 2 , psi 


PL f D 


= dimensionless pressure in a hydraulically fractured well, dimensionless 




= total compressibility, Lt 2 /m, psi" 1 




= average total compressibility, Lt 2 /m, psi" 1 . 
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37. The method of claim 36 wherein the straight line is defined by the equation: 
(y ap )n= b M + m M ( x a P )n , where 
_ (Pa)n ~ Par 



(y ap )n = 



[Pa(P)] jA -[Pa(P)]j 
(fa)j-(ta)jA 



, and 
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C op ]=C a j 

Cap2=Ca2 
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*Pa 
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(ta)n 

(%ap)n 

(yap)n 



3/2 



a first before-closure pressure-transient analysis adjusted variable, m/Lt 
psi 1/2 -cp 1/2 

a second before-closure pressure-transient analysis adjusted variable, 
m/Lt , psi -cp 

before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr, 

with adjusted pseudotime variable 

adjusted pressure variable difference, m/Lt 2 , psi 

adjusted reservoir pressure variable, m/Lt 2 , psi 

wellbore adjusted pressure variable, m/Lt 2 , psi 

time at timestep n, t, hr 

time at the end of an injection, t, hr 

adjusted time at timestep «, t, hr 

before-closure pressure-transient analysis adjusted variable, dimensionless 
before-closure pressure-transient analysis adjusted variable, dimensionless. 



38. The method of claim 37 wherein the first and second before-closure pressure- transient 
analysis variables are defined as: 



±- ;and 
<t>c t 
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wherein 

(j> = porosity, dimensionless 

B g = gas formation volume factor, dimensionless, bbl/Mscf 

B g = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L 2 t 2 , psi/ft 

w L = fracture lost width, L, ft. 




39. The method of claim 38 wherein the transforming step is iterated with a value of n varying 
from to a maximum value n max and for each couple of coordinates {(y ap )n, {*a P X} plot the 
graph (y ap )n versus (x ap )n to determine the slope m M , 
wherein 

ne = number of measurements that corresponds to the end of an injection 
IW = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure. 



40. The method of claim 39 wherein the permeability k is determined by: 

-a 

h (\A\ 9V9VO n9878V?zn I 



(141.2X2X0.02878)(24) 1 



5.615 



7 S f m M 



41. The method of claim 39 wherein the permeability k is determined by: 

-i2 



0)k = 



(141.2X2X0-02878X24) 1 



5.615 



7 S f m M 



wherein 

o) = natural fracture storativity ratio, dimensionless. 



42. The method of claim 34 wherein the injection fluid is a liquid, a gas or a combination 
thereof. 



43. The method of claim 42 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation 
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44. The method of claim 34 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 
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45. A method of estimating fracture-face resistance R 0 of porous rocks of a subterranean 
formation containing a compressible reservoir fluid comprising the steps of: 

(a) injecting an injection fluid into the subterranean formation at an injection pressure 
exceeding the subterranean formation fracture pressure; 

(b) shutting in the subterranean formation; 

(c) gathering pressure measurement data over time from the subterranean formation during 
shut-in; 

(d) transforming the pressure measurement data into corresponding adjusted pseudopressure 
data to minimize error associated with pressure-dependent reservoir fluid properties; and 

(e) determining the fracture-face resistance Ro of the subterranean formation from the adjusted 
pseudopressure data. 

46. The method of claim 45 wherein a plot of the adjusted pseudopressure data over time is a 
straight line with an intercept b/u a function of fracture-face resistance R 0 . 

47. The method of claim 46 wherein the adjusted pseudopressure data used in the transforming 
step are derived from the following equation: 

Ji c t pdp . 

(Pa )n = _ JLJL , wherein 

P Jo 





= average viscosity, m/Lt, cp 




= gas viscosity, m/Lt, cp 


p 


= pressure, m/Lt 2 , psi 


p 


= average pressure, m/Lt 2 , psi 


Pa 


= adjusted pseudopressure variable, m/Lt 2 , psi 


Pw 


= wellbore pressure, m/Lt 2 , psi 


PLfD 


= dimensionless pressure in a hydraulically fractured well, dimensionless 




= total compressibility, Lt 2 /m, psi" 1 




= average total compressibility, Lt 2 /m, psi" 1 . 



48. The method of claim 47 wherein the straight line is defined by the equation: 
(ya)n=™M(*a)n+t>M> where 
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/•„ \ _ (Pa)n-Par 
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wherein 

c al 
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Paw 

tn 

*ne 

(*a)n 

(y a )n 



= a first before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , 

^ -1/2^1/2 

psi -cp 

= a second before-closure pressure-transient analysis adjusted variable, 

2 n 2 J 12 -3/2 1/2 

m/Lt , psi cp 

= before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr 

= adjusted pressure variable difference, m/Lt 2 , psi 

= adjusted reservoir pressure variable, m/Lt 2 , psi 

= wellbore adjusted pressure variable, m/Lt 2 , psi 

= time at timestep «, t, hr 

= time at the end of an injection, t, hr 

= before-closure pressure-transient analysis adjusted variable, dimensionless 
= before-closure pressure-transient analysis adjusted variable, dimensionless. 



49. The method of claim 48 wherein the first and second before-closure pressure-transient 
analysis variables are defined as: 

c «' = i^ ;and 

5.615 _ Bg 
wherein 

(j> = porosity, dimensionless 
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B g = gas formation volume factor, dimensionless, bbl/Mscf 

B g = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L t , psi/ft 

w L = fracture lost width, L, ft. 

50. The method of claim 49 wherein the transforming step is iterated with a value of n varying 
from ne+l to a maximum value n max and for each couple of coordinates {(y fl )„, (x fl ) n }plot the 
graph (y a ) n versus (x a ) n to determine the intercept &a/, 

wherein 

ne = number of measurements that corresponds to the end of an injection 
n m ax = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure. 

51. The method of claim 50 wherein the fracture- face Ro is determined by: 

5.615 

*° = 141.2*(24) rpSftnebM 

52. The method of claim 45 wherein the injection fluid is a liquid, a gas or a combination 
thereof. 

53. The method of claim 52 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation. 

54. The method of claim 45 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 
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55. A method of estimating fracture- face resistance Ro of porous rocks of a subterranean 
formation containing a compressible reservoir fluid comprising the steps of: 

(a) injecting an injection fluid into the subterranean formation at an injection pressure 
exceeding the subterranean formation fracture pressure; 

(b) shutting in a zone of the subterranean formation; 

(c) gathering pressure measurement data over time from the subterranean formation during 
shut-in; 

(d) transforming the pressure measurement data into corresponding adjusted pseudopressure 
data and time into adjusted pseudotime data to minimize error associated with pressure- 
dependent reservoir fluid properties; and 

(e) determining the fracture-face resistance R 0 of the subterranean formation from the adjusted 
pseudopressure and adjusted pseudotime data. 

56. The method of claim 55 wherein a plot of the adjusted pseudopressure data over adjusted 
pseudotime data is a straight line with an intercept b M a function of fracture-face resistance R 0 . 

57. The method of claim 56 wherein the adjusted pseudotime and adjusted pseudopressure data 
used in the transforming step are respectively determined by: 





jj. = average viscosity, m/Lt, cp 

fi g = gas viscosity, m/Lt, cp 

p = pressure, m/Lt 2 , psi 

p = average pressure, m/Lt 2 , psi 

p a = adjusted pseudopressure variable, m/Lt 2 , psi 

p w = wellbore pressure, m/Lt 2 , psi 

PL f D = dimensionless pressure in a hydraulically fractured well, dimensionless 

c t = total compressibility, Lt 2 /m, psi' 1 

c t = average total compressibility, Lt 2 /m, psi" 1 . 
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58. The method of claim 57 wherein the straight line is defined by the equation: 
(y ap )n =b M +m M (x ap ) n , where 

\yap)n ~ r- r , 

\ a ap )n V" y'ne 



( d qp)j = 
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[Pa(P)] jA -[Pa(P)}j 
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wherein 

C a pi=C a j 
Cap2=Ca2 
dap 

Par 
Paw 

t n 

(ta)n 

(Xap)n 

(yap)n 



= a first before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , 
psi 1/2 cp 1/2 

= a second before-closure pressure-transient analysis adjusted variable, 

2/t 2.7/2 -3/2 1/2 

m/Lt , psi cp 

= before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr, 

with adjusted pseudotime variable 
= adjusted pressure variable difference, m/Lt 2 , psi 
= adjusted reservoir pressure variable, m/Lt 2 , psi 
= wellbore adjusted pressure variable, m/Lt 2 , psi 
= time at timestep n, t, hr 
= time at the end of an injection, t, hr 
= adjusted time at timestep «, t, hr 

= before-closure pressure-transient analysis adjusted variable, dimensionless 
= before-closure pressure-transient analysis adjusted variable, dimensionless. 



59. The method of claim 58 wherein the first and second before-closure pressure-transient 
analysis variables are defined as: 



c a\ s 



,zz~ j and 
</>c t 
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5.615 _ B g fjh 
wherein 

<t> = porosity, dimensionless 

£ g = gas formation volume factor, dimensionless, bbl/Mscf 

B g = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L 2 t 2 , psi/ft 

w L = fracture lost width, L, ft. 

60. The method of claim 59 wherein the transforming step is iterated with a value of n varying 
from ne+l to a maximum value n max and for each couple of coordinates {(y 0 p)n, (x ap ) n } plot the 
graph iy ap )n versus (x ap ) n to determine the intercept bM , 

wherein 

ne = number of measurements that corresponds to the end of an injection 
n ma x = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure. 

61. The method of claim 60 wherein the fracture- face R 0 '\s determined by: 

62. The method of claim 55 wherein the injection fluid is a liquid, a gas or a combination 
thereof. 

63. The method of claim 62 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation. 

64. The method of claim 55 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 
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65. A system for estimating physical parameters of porous rocks of a subterranean formation 
containing a compressible reservoir fluid comprising: 

(a) a pump for injecting an injection fluid into the subterranean formation at an injection 
pressure exceeding the subterranean formation fracture pressure; 

(b) means for gathering pressure measurement data from the subterranean formation during a 
shut-in period; 

(c) means for transforming the pressure measurement data into adjusted pseudopressure data to 
minimize error associated with pressure-dependent reservoir fluid properties; and 

(d) means for determining the physical parameters of the subterranean formation from the 
adjusted pseudopressure data. 

66. The system of claim 65 wherein the determining means comprises graphics means for 
plotting a graph of the adjusted pseudopressure data over time, the graph being a straight line 
with a slope % and an intercept bM wherein yyim is a function of permeability k and b^ is a 
function of fracture-face resistance Ro* 

67. The system of claim 66 wherein the adjusted pseudopressure data is defined by the 
following equation: 



P 



^ , wherein 



g t 



{i = average viscosity, m/Lt, cp 

jj g = gas viscosity, m/Lt, cp 

p = pressure, m/Lt 2 , psi 

p = average pressure, m/Lt 2 , psi 

p a = adjusted pseudopressure variable, m/Lt 2 , psi 

p w = wellbore pressure, m/Lt 2 , psi 

PL f D = dimensionless pressure in a hydraulically fractured well, dimensionless 

c t = total compressibility, Lt 2 /m, psi" 1 

c t = average total compressibility, Lt 2 /m, psi' 1 . 

68. The system of claim 67 wherein the straight line is defined by the equation: 
(y a )n = m M ( x a)n +&M > where 
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/„ x _ (Pa)n-Par 
\sa in - ~~~~ I — I — 



[Pa(P)]j-] -[Pa(P)]j 



, and 



(*a)n s 



( d a)ne+2 



V2 



(d a ) n 



+ z 

j=ne+3 
c a2 



\ a a)n l ne 



1- 



^1 t m¥\ ^ 
*n J 



wherein 

c a\ 

c aZ 

d a 

&Pa 

Par 

Paw 

'« 

*ne 

(y a ) n 



a first before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , 

•1/2 1/2 

psi -cp 

a second before-closure pressure-transient analysis adjusted variable, 

2 /T 2.7/2 -3/2 1/2 

m/Lt , psi cp 

before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr 

adjusted pressure variable difference, m/Lt 2 , psi 

adjusted reservoir pressure variable, m/Lt 2 , psi 

wellbore adjusted pressure variable, m/Lt 2 , psi 

time at timestep n, t, hr 

time at the end of an injection, t, hr 

before-closure pressure-transient analysis adjusted variable, dimensionless 
before-closure pressure-transient analysis adjusted variable, dimensionless. 



69. The system of claim 68 wherein the first and second before-closure pressure-transient 
analysis variables are defined as by the equations: 



and 



5.615 c B g Uig 
wherein 

^ = porosity, dimensionless 
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B g = gas formation volume factor, dimensionless, bbl/Mscf 

B g = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L 2 t 2 , psi/ft 

w L = fracture lost width, L, ft. 

70. The system of claim 69 wherein the transforming means iterates the transformation of each 
adjusted pseudodata with a value of n varying from ne+l to a maximum value n max , and 
wherein the graphics means plots the graph (y a ) n versus (x 0 ) n to determine the slope m M and the 
intercept b M , 

wherein 

ne = number of measurements that corresponds to the end of an injection 
n max = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure. 

71. The system of claim 70 wherein the permeability k and the fracture-face R 0 are determined 
by the following equations: 

; and 



k = 



(141.2)(2)(0.02878)(24) 1 



5.615 r p Sftn M 



*° = 141.2.(24) r P S f f ne b M 



72. The system of claim 70 wherein the permeability k and the fracture-face ^are determined 
by the following equations: 

; and 



cok = 



(141.2)(2)(0.02878)(24) 1 



5.615 r p Sfm M 



5 - 615 

^ = 141.2*<24) r P S f 
wherein 

o) = natural fracture storativity ratio, dimensionless. 

73. The system of claim 65 wherein the injection fluid a liquid, a gas or a combination thereof. 
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74. The system of claim 73 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation. 

75. The system of claim 65 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 
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76. A system of estimating physical parameters of porous rocks of a subterranean formation 
containing a compressible reservoir fluid comprising: 

(a) a pump for injecting an injection fluid into the subterranean formation at an injection 
pressure exceeding the subterranean formation fracture pressure; 

(b) means for gathering pressure measurement data from the subterranean formation during a 
shut-in period; 

(c) means for transforming the pressure measurement data into adjusted pseudopressure data 
and time into adjusted pseudotime data to minimize error associated with pressure-dependent 
reservoir fluid properties; and 

(d) means for detecting characteristics of the evolution in the adjusted pseudopressure data over 
adjusted pseudotime data to determine the physical parameters of the subterranean formation. 

77. The system of claim 76 wherein the detecting means comprises graphics means for plotting 
the evolution of the adjusted pseudopressure data over adjusted pseudotime data, the evolution 
being a straight line with a slope m M a function of permeability k and an intercept b M a function 
of fracture-face resistance Ro. 

78. The system of claim 77 wherein adjusted pseudotime and adjusted pseudopressure data are 
respectively determined by the equations: 





p. = average viscosity, m/Lt, cp 

/j g = gas viscosity, m/Lt, cp 

p = pressure, m/Lt 2 , psi 

p = average pressure, m/Lt 2 , psi 

p a = adjusted pseudopressure variable, m/Lt 2 , psi 

p w = wellbore pressure, m/Lt 2 , psi 

PL f D = dimensionless pressure in a hydraulically fractured well, dimensionless 

c t = total compressibility, Lt 2 /m, psi' 1 
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= average total compressibility, Lt 2 /m, psi" 1 . 



79. The system of claim 78 wherein the straight line is defined by the equation: 

(yap )n= b M + m M (*ap )n > where 

/„ x _ (Pa)n ~ Par 
\sap )n ~ 



v ap>n — ( . J — J , 



[Pa(P)]j-rlPa(P)]j 
(ta)j-(QjA 



, and 



( x ap)n - 



c ap\ 



(dap)ne+2 



( d ap)n 



V2 



t n t } 



ne 



+ 



j=ne+3 



c ap2^(^n 



id ) t^ 2 



( d ap)n 



(fa)n 



wherein 

C a p]=C a j 

dap 

*Pa 

Par 

Paw 

*n 

*ne 

Oa)n 

(Xap)n 

(yop)n 



a first before-closure pressure-transient analysis adjusted variable, m/Lt 3/2 , 

•1/2 \I2 

psi *cp 

a second before-closure pressure-transient analysis adjusted variable, 

2/t 2/7/2 -3/2 1/2 

m/Lt , psi xp 

before-closure pressure-transient analysis adjusted variable, m/Lt 3 , psi/hr, 

with adjusted pseudotime variable 

adjusted pressure variable difference, m/Lt 2 , psi 

adjusted reservoir pressure variable, m/Lt 2 , psi 

wellbore adjusted pressure variable, m/Lt 2 , psi 

time at timestep n, t, hr 

time at the end of an injection, t, hr 

adjusted time at timestep n, t, hr 

before-closure pressure-transient analysis adjusted variable, dimensionless 
before-closure pressure-transient analysis adjusted variable, dimensionless. 



80. The system of claim 79 wherein the first and second before-closure pressure-transient 
analysis variables are defined as: 
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t 



c a \ = x hr& ;and 

5.615 J*g\h . 
24 7 (B g ) ne \<t>c t 

wherein 

^ = porosity, dimensionless 

£ g = gas formation volume factor, dimensionless, bbl/Mscf 

2? g = average gas formation volume factor, dimensionless, bbl/Mscf 

Sf = fracture stiffness, m/L 2 t 2 , psi/ft 

w L = fracture lost width, L, ft. 

81. The system of claim 80 wherein the transforming means iterates the transformation of each 
adjusted pseudodata with a value of n varying from ne+l to a maximum value n max , and 
wherein the graphics means plots the graph (y a ) n versus (x a ) n to determine the slope m M and the 
intercept b M , 

wherein 

ne = number of measurements that corresponds to the end of an injection 

n m ax = corresponds to the data point recorded at fracture closure or the last recorded 
data point before induced fracture closure. 

82. The system of claim 80 wherein the permeability k and the fracture-face R 0 are determined 
by the equations: 

i2 



k = 



(141.2)(2)(0.02878)(24) 1 



5.615 



'P 



S f m M 



; and 



^T^^) r p s f tnebM 



83. The system of claim 80 wherein the permeability k and the fracture-face R 0 are determined 
by the equations: 

-|2 

6)k = 



(141.2)(2)(0.02878)(24) 1 



5.615 r p Sfm M 
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Ro = m2 5 m r r S f' nebM '> 

wherein 

o) = natural fracture storativity ratio, dimensionless. 

84. The system of claim 76 wherein the injection fluid is of a liquid, a gas or a combination 
thereof 

85. The system of claim 84 wherein the injection fluid contains desirable additives for 
compatibility with the subterranean formation. 

86. The system of claim 76 wherein the reservoir fluid is a liquid, a gas or a combination 
thereof. 
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